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Abstract 

The reduction of the n per rev. pitch-, roll- and vertical 
vibrations of an n-bladed rotor by n per rev. sinusoidal variations 
of the collective and cyclic controls is investigated. The 
numerical results presented refer to a four-bladed, 7.5-foot model 
and are based on frequency response tests conducted under an 
Army-sponsored research program. The following subjects are 
treated: 

• Extraction of the rotor transfer functions (.073R hub 
flapping and model thrust versus servo valve command, 
amplitude and phase) 

• Calculation of servo commands (volts) required to 
compensate .073R hub flapping (3P and 5P) and 
model thrust (4P) 

• Evaluation of the effect of the vibratory control inputs 
on blade loads 

• Theoretical prediction of the root flapbending 
moments generated by o to 5P perturbations of the 
feathering angle and rotor angle of attack. 

Five operating conditions are investigated covering advance 
ratios from approximately 0.2 to 0.85. The feasibility of vibra- 
tion reduction by periodic variation on conventional controls is 
evaluated. 


For several operating conditions covering advance ratios 
from approximately 0.2 to 0.85, the control inputs required to 
counteract the existing 4P pitch, roll and vertical vibrations are 
calculated. The investigations are based on experimental vibra- 
tion and response data. As the tests were part of and added on to 
a larger hingeless rotor research program, only a few operating 
conditions with essentially zero tip path plane tilt were investi- 
gated because of limited tunnel time. At the test rotor speed (500 
rpm) the rotor blade mode frequencies were 1.34P, first flapping, 
6.3P, second flapping, and 3.6P, first inplane. 


This work was conducted under the sponsorship of the 
Ames Directorate of the U. S. Army Air Mobility R&D Lab- 
oratory under Contract NAS2-7245. The authors gratefully 
acknowledge the assistance of Mr. David Sharpe, the AMRDL 
Project Engineer, and Messrs. R. London and G. Watts of 
Lockheed in conducting the experimental portion of this work. 


It should be noted that there was no instrumentation to 
measure the vibratory pitching and rolling moments. These 
moments were obtained by properly adding up the flap-bending 
moments of the four blades at 3.3 in. (0.073R) which were 
measured separately. This means, the effects of the inplane 
forces, vertical shear forces and blade torsion have been ignored. 
These are important influences in current hingeless rotor designs. 
The inplane 3P and !>P shear forces are of particular 
interest. However, the experimental data obtained for a model 
hingeless rotor system provides the beginning of at least a partial 
data base for the investigation of vibration attenuation of such 
systems through periodic variation of conventional controls. 

Generally speaking, the control inputs required for flapping 
(hub moment) sourced vibration elimination are smaller or about 
of the same magnitude as those used for the frequency response 
tests. Their amplitudes lie, depending on flight condition and 
advance ratio, between 0.2 and 3 degrees. With the exception of 
the /i = 0.851 case, for which the results are somewhat in doubt 
(the response tests to lateral cyclic pitch and the corresponding 
baseline data were inadvertently run with 0.3-degree collective 
pitch differential), the control inputs required for vibration re- 
duction drastically reduce the 3 and 5P, and have only a minor 
effect on the 2P flexure flap-bending moments. Chord-bending 
moments and blade torsion generally increase. 

The theoretical predictions mentioned refer to forced- 
response influence coefficients. They are based on the first two 
flapping modes. The blade root flap-bending moments (OP 
through 5P) which result from unit perturbations of blade 
feathering angle and rotor angle of attack have been calculated. 
The solution provides for intermode coupling through the 17th 
harmonic by analytic solution of the two-degree-of-freedom 
system, utilizing constant coefficient and loading descriptions 
over ten-degree azimuth sectors. In each solution case, the rotor 
reached steady-state motion in eight revolutions. In that time the 
least converging second mode flapping motion converged to a 
minimum of four significant figures. 

Evaluation of the test data reveals two types of short- 
comings, which should be avoided in future tests. First, the data 
given are based on a single test and have not been verified. 

Second, in some cases, the baseline and frequency response tests 
were not run successively. 

From the data available, the approach is promising, 
especially for the low and medium advance ratio range. At higher 
advance ratios (p~ 0.8), the control inputs required for vibration 
reduction may become prohibitive. 

Notation 


A, B 


Presented at the AHS/NASA-Ames Specialists* Meeting on 
Rotorcraft Dynamics, February 13-15, 1974. 
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quantities describing cos 4^ and sin 4 ^components 
of actuator input for frequency response tests, volt, 
see Table II and Equation (1) 


C, D quantities describing responses to A and B, in.-lb 

and lb, respectively, see Equation (1) 

E, F, G, H blade loads due to unit actuator input, in.-lb/ volt, 
see Equation (13) 

Kj . . . Kjg gains of rotor response, see Table I 

m calculated flapbending moment at 3.3 in., in.-lb, 

m = m Q + Xm ns sin n^ + Sm nc cos n\fr 

M, L, T 4P vibratory pitching moments, rolling moments 
and thrust variations, in.-lb and lb, respectively; 
subscript e denotes existing vibrations to be com- 
pensated, subscript control describes effects of 
oscillatory control inputs. 

M e = M s sin 4^ + M c cos 4^ 

L e = L s sin 4 ^ + -L C cos 4^ 

T e = T s sin 4\l' + T c cos4* 

^nominal nominal collective pitch, degrees 

0 Q , 0 S , 0 C oscillator inputs for collective, longitudinal 
and lateral cyclic pitch, volt 


the vibratory blade loads. These studies are based on frequency 
response tests conducted on a 7.5 foot-diameter, four-b laded, 
hingeless rotor model, the results of which are published in 
Appendixes C and D of Reference 1. The subject matter covered, 
apart from the items listed below, is an abridged version of these 
appendixes. 

Other subjects treated are (a) the calculation of blade loads, 
based on test data, due to vibratory control command applica- 
tions; (b) the theoretically determined eigenvalues, at 10-degree 
azimuth intervals, of the first and second flapping modes, at 
v = 0.191, 0.45 and 0.851 ; (c) the computed single-blade root 
flap-bending moment, Sta 0, harmonic influence coefficients 
at At = 0. 1 9 1 , 0.45 and 0.851; and (d) a limited comparison of the 
theoretical loads with experiments. 

The general case of vibration control will include the effects 
of lateral and fore-and-aft shear forces at blade passage frequency. 
These forces can be as influential as the pitch and roll moment 
and thrust oscillations in causing fuselage vibrations. Thus, in 
general, five rotor vibratory inputs are to be controlled by mani- 
pulation of three controls. Although the five vibratory inputs 
cannot be nulled individually with three controls, their combined 
contribution to the fuselage vibration can be controlled. Thus, 
the general application will involve control of fuselage vibration 
at three points; say two vertical vibrations and one roll angular 
vibration. This general application implies the use of adaptive 
feedback controls. Although the present paper is limited to the 
more simple case outlined herein, the general application to the 
control of any three suitable quantities will be apparent. 


6 o = e os sin 4 ^ + e oc cos 4 ^ 

0 S = ^ss s * n 4 ^ + ^sc cos 4 ^ 

0 C = 0 CS Sin 4^ + 0 CC cos 4^ 

T 1 ’ • • r 18 la 8 angles of response, degrees, see Table I 


Although prior investigations of the use of higher harmonic 
pitch control on teetering and offset hinge rotors have been con- 
ducted to investigate improved system performance and also for 
vibration attentuation (References 2 , 3 and 4 ), this is believed 
to be the first experimental and theoretical hingeless rotor study 
of the use of periodic variation of conventional controls for 
vibration attentuation. The use of 2P feathering to improve rotor 
performance is not included as part of this work. 


a 


rotor angular velocity, sec"* 


Transfer Functions Involved 


0 


azimuth position of master blade, rad 


r 

^RM = 
a o 


Blade Root Moment, STA(o) 
ttR 3 p(£2R)^ao 


where 


a = 5.73 


As a distinction must be made between control applications 
in phase with sin 4^ and cos 4^, there are six control quantities 
available, i.e., 0 QS , 0 QC , 0 sg , 0 SC , 0 C§ and 0 CC , to monitor the 
pitching moments, rolling moments and vertical forces. This 
means the dynamic system investigated, which consists of rotor, 
control mechanism and oscillators used, is characterized by 18 
gains K p and lag angles r p . The subscripts p (p = 1 through 18) 
are defined by Table I. 


P = 0.002378 slugs/ft 3 

o = 0.127 

“Compensating Control Inputs” define those which reduce 
the existing 4P pitching moments, rolling moments and vertical 
forces of a given flight condition to zero. 

The analysis deals with the concept of vibration reduction 
by oscillatory collective and cyclic control applications. Several 
related aspects of this problem are treated. The foremost are the 
determination of the proper control inputs and their effect on 


TABLE I 

GAINS AND LAG ANGLES OF RESPONSE 
TO OSCILLATORY CONTROL APPLICATIONS 



*os 

6 oc 

0 ss 

®sc 

e cs 

e cc 

M 

K l r l 

k 2 t 2 

K 3 t 3 

K 4 r 4 

K 5 t 5 

K 6 '6 

L 

K 7 r 7 

K 8 t 8 

k 9 t 9 

K 10 t 10 

K irn 

K 12 t 12 

T 

K 13 t 13 

K 14 t 14 

K 15 t I5 

K 1 6 r 1 6 

K 1 7 r l 7 

K 18 r 18 
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As indicated, K 3 is defined as the amplitude ratio M 10^ and 
r 3 is the lag angle of M with respect to 0^. For convenience, the 
dimensions used are identical with those of the computer output, 
i.e., oscillator voltage for input, in.-lb for M and L, lb for the 
thrust variation T. This means the dimensions of Kp are 

Kj through Kj 2 in.-lb/volt 

K j 3 through K j g lb/volt 


See also Figure 1 which shows the oscillatory pitching moments 
due to combined ^ss and 6 SC control applications. The moments 
generated are presented by rotating vectors where cos 4 ^ is posi- 
tive to the right and sin 4 $ positive down. This means, the vector 
positions shown refer to $ = 0. By definition, the quantities Ry 
characterize the responses in phase with the excitation and Ij : 
those out of phase. The latter are positive if the response leads. 
As indicated, there are altogether four responses involved which 
are combined to the resultant M. 


The phase angles r p are given in degrees, r p is positive if the re- 
sponse lags. 

Although the investigations deal exclusively with 4P control 
variations, some general remarks may be in order. The general 
case involves sinusoidal collective and cyclic control variations 
with the frequency nn where n can be any positive number. 

If n is an integer, the rotor excitations repeat themselves 
after each rotor revolution which means that the responses of 
each revolution are identical. This is true for any number of rotor 
blades but does not necessarily mean that all blades execute 
identical flapping motions. The latter is true only if n equals the 
number of rotor blades or is a multiple of the blade number. 

Only for these cases does a truly time independent response with 
invariable amplitude ratios K and lag angles r exist. 

Extraction of Gains and Lag 
Angles from Experiments 


As for all response tests conducted, the oscillator input con- 
tained both sin 4^ and cos 4^-components; always two amplitude 
ratios K and two lag angles r are involved. Therefore, each time a 
set of two tests must be evaluated. According to Table II, the 
input is characterized by the quantities Aj Bj A 2 B 2 and the 
response by Cj Dj C 2 D 2 . 

If the rotor responds to cos 4^ excitations with the gain Kj 
and the lag angle Tj (j = even number) and to sin 4 ^ excitations 
with Kj and Tj (i = odd number), input and output are related by 
the equations 


Aj Kj cos (4^- 
A 2 Kj cos (4tf - 


7j) + Bj Kj sin (4^ - Tj) - Cj cos 4 $ 

+ Dj sin 4^ 

Tj) + B 2 Kj sin (4^ - 7j) = C 2 cos 4^ 

+ D 2 sin 4^ J 


( 1 ) 



Figure 1 . Vector Diagram Showing Pitching Moment 
Due to 0 ^ and 0 SC Control Applications 

Inserting Equation (2) into Equation (1) leads to 


a 1 D 2 ‘ A 2 D 1 


Aj B 2 ■ A 2 B j 

j A 1 C 2 ' A 2 C 1 
* A 1 ®2 * A 2®1 


(3) 


TABLE II 

INPUT AND OUTPUT NOTATIONS 


Test 

Input 

Response 

#1 

Aj cos 4^ + Bj sin 4^ 

Cj cos 4^ + Dj sin 4^ 

#2 

A 2 cos 4 ^ + B 2 sin 4^ 

C 2 cos 4^ + D 2 sin 4^ 


tan?j= Ilj/R|l OCTjOr /2 J 


and 

c 1 b 2 - b 1 c 2 
*3 A j B 2 - A 2 Bj 


To calculate the unknowns Kj Kj Tj and Tj, a component 
analysis is used. The gains Kj Kj are expressed as 

Kj = (R 2 + Ip 1/2 

Kj = (R 2 + I 2 ) 1/2 


^ ( 2 ) 


BjD 2 * B 2 D j 
1 a 1 b 2 ' A 2 B 1 

tan Tj = I Ij/Rj I 0 < Tj <tt/2 


( 4 ) 
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In both cases 


r = +T for 

R> 0 

I<0 

= - 7 

R > 0 

I>0 

= 7 r + r 

R< 0 

I>0 

= 7 r - r 

R< 0 

ko 


Check of Calculated Kj Kj 7j and rj Values 


If so desired, Equation ( 1 ) can be used to check the calcu- 
lated values of Kj Kj r- and r • Splitting up these equations into 
sin 4 ^ and cos 4 y components leads to the following four 
expressions which must be satisfied 


Aj Kj cos rj - Bj Kj sin r* = Cj 
Aj Kj sin 7j H-BjlCi COS 7 j = Dj 
A 2 Kj cos 7j - B 2 Kj sin 7j = C 2 
A 2 Kj sin tj r B 2 Kj cos t* — D 2 


( 5 ) 


Oscillatory Control Inputs Required 


The six oscillator inputs available have to be selected so that 
their responses satisfy the requirements, whatever they may be. 
By definition, the vibratory control inputs result in the following 
pitching moments, rolling moments and vertical forces (n = 4): 

M control = + 0 os K l sin(n^- r,) 

+ d QC K 2 cos (n0 - rj) 

+ 0 SS K3 sin (nP - T3) 

+ 0 SC K 4 cos (n^- t 4 ) 

+ 0 CS Kg sin (ntp - »g) 

+ 0 C c K 6 cos ( n ^- T 6> (6) 


^control + ^os ^7 s ' n ' T 7^ 

+ 9 oc Kg cos(n ip- r 8 ) 

+ 8 SS Kg sin (np- Tg) 

+ 0 sc K, o cos (n<P- r 1Q ) 

+ 0 CS Kj j sin (nip- t u ) 

+ 0 cc K 12 cos(nl ^' t 12> ( 7 ) 


^control + ^os ^13 s * n * n ^ / ' T 13 ) 

+ 0 QC K 14 cos (nil'- r, 4 ) 

+ 6 SS K 15 sin (nil'- r,g) 

+ e sc K 16 C0S < n ^- t 16 > 

+ 0 cs K 17 sin (n<P- r 1? ) 

+ 0 C c K i8 co s(n'/'- ri 8 ) (8) 


M C ontrol = -M s sin " M c cos 

^control “ -L s s * n ^ ' ^c cos ^ 

T J. . » = -T sin All/ - T A\L/ 

V/Uiitroi * s A C — “ ■ ' 


( 9 ) 


To reduce the existing vibrations, the moments and forces 
generated must counteract M e , L e and T e , i.e., 

Equations 6 through 9 lead to six linear equations, (10), 
for the unknowns d os , 0 OC , 0 SS> 6 SC , 6 CS and 8 CC 

Effect on Blade Loads 


An objective of the investigations is to determine the effect 
of the compensating control input on the blade loads, i.e., on the 
following measured quantities: 

• flapbending at 3.3 in. 

• flapbending at 1 3. 1 5 in. 

• chordbending at 2.4 in. 

• torsion at 9.28 in. 

In all cases the 2 to 5 P content of the loads is of interest. 
The first task is to determine from the response tests the contri- 
bution of each of the six possible 4 P control inputs to these 
loads. Again, two sets of data are required. The vibratory control 
applications used and the resulting n^ 1 harmonic of the load con- 
sidered are written as follows: 


Test 

Input 

Resulting Load (in.-lb) 

#1 

Aj cos 4 ^ + Bj sin 4^ 

C n j cos n^ + D n j sin n^ 

#2 

A 2 cos 4 ^ + B 2 sin 4^ 

C n 2 cos n^ + D n2 sin n'/' 



*+Kj COS Tj 

+K 2 sin r 2 

+K3 cos 7^ 

+K4 sin r 4 

-Kjsin 7j 

+K 2 cos r 2 

-K3 sin 

+K4 cos 74 

+ Ky COS T'j 

+Kg sin rg 

+K9 COS Tg 

+K j 0 sin r 1Q 

-Kj sin Tj 

+Kg COS Tg 

-Kg sin Tg 

+ K 1 0 cos T 10 

+Kj 3 cos t ,3 

+K j4 sin Tj4 

+K15 cos t , 5 

+K 16 sm r, 6 

_- R 1 3 sin r 13 

+Kj 4 cos 7 j 4 

-Kigsin T,g 

+Ki6 c °s r, 6 


+K5 COS 73 

+K 6 sin r 6 


^os 



-K5 sin 73 

+K 6 cos t 6 


^oc 


-M c 

+Kj j cos 7 j j 

+K j 2 sin t j 2 


^ss 


-h 

-Kn sin rj j 

+ K 1 2 cos T 1 2 


^sc 


L c 

+K17 cos 7 j 2 

+ K lg sin r lg 


&cs 


-T s 

-Kj 2 sin 7 j 2 

+ K ] 8 cos t 18 


0 

1 


T c 


( 10 ) 
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If nonlinear effects are excluded, the n per rev load vari- 
ation due to unit control application in phase with 

(a) cos 4^ amounts to (E n cos n'l' + F n sin n'P ) 

(b) sin 4* (G n cos n'/' + H„ sin iWO 

In these expressions 


► 


( 12 ) 


= 


F„ = 


- 


»n = 


®2 C nl 

B lC n 2 

a 1 b 2 " 

a 2 B 1 

B 2 D nl 

‘ B l D n2 

A 1 B 2 ” 

a 2 B 1 

A l C n2 

* A 2 c nl 

a 1 b 2 ‘ 

a 2 B 1 

A lPn2 

- a 2^i1 

a 1 b 2 ' 

a 2 B 1 


H3I 


used. 

are 


If 0t s , 0t c (£ = o, s, c) denote the vibratory control inputs 
, the increments of the n^ 1 harmonic of the load considered 


(Aload) n = (0g c E n + 0^ G n ) cos n* 

+ ( %c F n + 0 £s H n )sin nJ/ 
Evaluation of Experiments 


(14) 


Flight Conditions Investigated 


The methods outlined in the previous sections are applied to 
the following five operating conditions for which test data are 
available: 


TABLE III 

OPERATING CONDITIONS INVESTIGATED 


M 

^nominal 

a 

Cj/o 

0.191 

12° 

-5° 

0.102 

0.239 

4 

-5 

0.028 

0.443 

4 

-5 

0.011 

0.849 

10 

-5 

-0.005 

0.851 

4 

-5 

-0.013 


In all cases the shaft angle of attack is ft = -5° and the rotor is 
trimmed so that essentially aj = b| =0. As can be seen, the tests 
cover the advance ratio range from approximately M = 0.2 to 
M= 0.85. The case 0.191 is characterized by 0 nominal = 
and C T /a = 0.102, the latter figure indicates a relatively high 
specific loading. In contrast, at the advance ratios M = 0.849 
and 0.851 the rotor is practically unloaded, i.e., no steady lift- 
ing force is generated. The 4P vibrations associated with the 
various test conditions are listed in Table IV. The moments are 
given in inch-pounds and the vibratory forces in pounds. 


These moments were obtained by properly adding up the 
flap-bending moments of the four blades at 3.3 in. which were 
measured separately. This means, the effects of the in-plane 
forces, vertical shear forces and blade torsion have been ignored. 


TABLE IV 

VIBRATORY MOMENTS AND FORCES 
TO BE COMPENSATED 


M 

0.191 

0.239 

0.443 

0.849 

0.851 

Ms 

0.3805 

-1.7207 

2.6149 

20.0483 

3.5349 

M c 

- 0.5301 

-0.4113 

- 0.5208 

-4.5724 

-8.4341 

h 

12.2080 

1.3725 

- 6.7626 

9.4647 

-10.5154 

L c 

2.2180 

-1.9145 

- 3.7399 

-31.1214 

-17.2626 

T s 

0.1979 

-0.1089 

0.0304 

1.9247 

0.8838 

T c 

-0.2013 

- 0.0865 

0.0556 

- 0.0048 

- 0.8626 


Gains and Lag Angles 


The rotor response characteristics are calculated by applying 
equations (2, 3, 4) to the test data available. The results 
available are listed in Table V. As pointed out previously, the 
values given include the effect of the actuator used. Some 
general statements can be made. It is obvious that for n = 0. 
the gain and lag angle of the responses to sin 4^ and cos 
4 v -type control applications must be the same. For ^ ^ 0 this is 
no longer true, and one would expect that the spread between 
KjKj and t-tj (see equations (3), (4)) widens with increasing 
advance ratio. Further, according to classical rotor theory which 
neglects blade stall, the nominal collective pitch setting has no 
effect on the frequency response characteristics. 

Generally speaking, the KjKj and tjtj values given in 
Table V differ very little. It appears however, that at higher 
advance ratios (compare columns for m = 0.849 and 0.851) the 
collective pitch has a larger effect than anticipated. It is also 
possible that the error of the baseline data described in the sum- 
mary may play a role. 

Oscillator Inputs Required 

Equation (10) is used to calculate the inputs required to 

(a) generate unit amplitudes of pure pitching moments, 
rolling moments and vertical forces and 

(b) compensate the existing vibrations 

The results are given in Tables VI and VII. They show that, as 
to be expected, the oscillatory inputs required for vibration 
reductions generally increase with increasing advance ratio. 
Surprisingly, the rotor collective pitch setting seems to play a 
larger role than the steady lift generated. See also Table VIII 
which summarizes the results obtained and lists the operating 
conditions investigated in the order of decreasing vibrations. 

The first column shows the relative magnitude of the vibratory 
moments generated and the last column the approximate 
amplitude of the blade pitch variation required to compensate 
the vibrations. The amplitude of the pitch variation produced 
per volt oscillator input changes with the control loads and 
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TABLE V 

GAINS AND LAG ANGLES DERIVED FROM EXPERIMENTS 
(Kp - in.-lb/volt, T p - degrees) 



P = 0. 

191 

M = 0. 

239 

n = 0. 

443 

M= 0.849 

M=0. 

851 

p 

s 

T P 

S 

t p 

Kp 

t p 

K p 

t p 

K p 

t p 

1 

5.617 

42.3 

1.099 

125.6 

2.236 

120.5 

4.798 

72.0 

4.094 

116.5 , 

2 

6.126 

44.0 

1.141 

149.1 

2.791 

129.3 

4.787 

72.6 

3.487 

135.6 

3 

17.571 

-9.6 

52.416 

-30.1 

42.237 

-28.7 

18.537 

- 19.8 

43.319 

-5.1 

4 

26.019 

-45.4 

47.991 

-37.3 

40.073 

-30.1 

20.329 

-41.5 

37.081 

12.7 

5 

30.696 

155.7 

59.416 

182.9 

45.186 

188.4 

33.002 

183.4 

26.170 

214.2 

6 

32.505 

181.7 

77.408 

193.2 

61.144 

180.8 

21.085 

180.0 

38.661 

184.5 

7 

2.856 

136.0 

4.246 

81.9 

8.166 

86.5 

2.472 

102.1 

10.097 

93.1 

8 

i .507 

98.4 

5.083 

67.1 

8.077 

66.9 

3.412 

144.7 

7.979 

62.9 

9 

35.384 

213.4 

59.420 

198.8 

43.846 

181.4 

44.506 

200.5 

48.081 

176.2 

10 

41.674 

185.8 

51.280 

198.6 

39.383 

195 7 

48.473 

201.0 

40.850 

187.7 

11 

45.953 

116.6 

76.875 

108.3 

78.512 

101.8 

67.268 

134.4 

88.540 

94.7 

12 

61.589 

131.5 

86.361 

99.3 

80.995 

95.7 

61.288 

141.5 

90.934 

95.3 

13 

6.879 

45.6 

5.420 

51.4 

8.928 

39.2 

8.188 

35.8 

9.340 

38.5 

14 

7.211 

43.7 

6.195 

46.4 

8.999 

35.9 

8.906 

36.1 

; 9.651 

35.6 

15 

6.635 

245.2 

4.275 

205.9 

2.571 

195.2 

5.976 

215.0 

3.623 

184.0 

16 

6.033 

218.3 

3.962 

208.1 

3.123 

188.7 

4.775 

229.5 

1.977 

185.4 

17 

13.000 

127.3 

7.596 

94.3 

7.632 

76-7 

13.261 

133.1 

11.188 

86.9 

18 

10.057 

128.6 

8.176 

97.4 

8.381 

92.2 

7.953 

126.3 

11.101 

90.7 


the type of control (0 O ,0 s , 0 C ) used. Therefore, the con- 
version factor varies and the last column of Table VIII is 
given only to indicate the approximate amplitudes involved. 

With one exception, the vibratory control applications re- 
quired were smaller than those used for the frequency response 
tests. The exception is the case with the highest vibration level 
encountered for which the compensating controls required were 
approximately 15 to 20% higher than the inputs used for the 4P 
frequency response tests. 

Blade Loads 


The calculation of the effect of the compensating control 
inputs on the blade loads is based on Equations (13) and ( 1 4). The 
first step is to calculate, for each specific case, the quantities E n 
through H n (n = 2, 3,4, 5). See Table IX which refers to /i =0,849 
and lists the sin mjjand cos n+ components of the various 
loads due to unit control (volt) application. The table shows, for 
instance, that at the advance ratio m = 0.849, a ±1 volt variation 
of 0 SS produces 3P chordwise bending moments of the magnitude 

(-91.77 sin 3^ +7. 1 5 cos 3$) in.-lb 


As the control inputs required for vibration reduction have been 
previously calculated, their effects on the blade loads can be de- 
termined by adding up the various contributions. The reader is 
referred to Table X which applies to the flapbending moment 
at 3.3 in. for the case m= 0.849. Given are the original loads 
without vibratory control application, the individual contribu- 
tions and the sum. The last column shows the amplitudes with- 
out and with compensating control input. A summary of the 


loads is represented in Table XL Generally speaking, chord- 
bending, blade torsion and the 4P flap-bending moments of the 
root flexure increase with increasing advance ratio. The 3 and 5P 
flap-bending moments of the flexure are, by nature, reduced and 
the 2P flap-bending moments are least affected. From the limited 
data available, it appears that the 4P chordwise- and 5P torsion 
moments may be the critical load for this configuration, 
inasmuch as the natural frequencies are close to these 
values. 


As mentioned previously, it is assumed here that 
the pitching and rolling moments are solely caused by the flap- 
bending moments of the root flexure which were individually 
measured and properly combined by a sin-cos potentiometer. 

This means, the only source for the troublesome 4P moments in 
the nonrotating system are the 3 and 5P flap-bending moments at 
3.3 in. For four identical blades, it follows that elimination of 
the 4P pitching and rolling moments requires that the sin 3^, cos 
3^, sin 5^ and cos 5^ components of the flap-bending moments 
at 3.3 in. are reduced to zero. As the four blades behave dif- 
ferently, this ideal condition will practically never be fulfilled. 

In the preceding paragraphs the flapbending moment of a 
specific blade, with consideration of the compensating control 
input, was calculated. To a certain extent, these predicted loads 
can be used as an independent check. As an example, the case 
M- 0.849 is treated. According to Table IV the amplitudes of 
the 4P pitching and rolling moments to be compensated are 


M = 20.56 in.-lb 
L = 32.52 in.-lb 


(15) 
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The calculated 3 and 5P flap-bending moments with considera- The amplitudes of the resulting 4P pitching and rolling moments 
tion of the compensating control input amount to (see Table VII), are 


m 3s = 0.6233 in.-lb 
m 3c - -1.1833 
m^ s — -1.9266 
m 5c = 0.3099 


(16) 


M = 3.14 in.-lb 
L = 5.91 in.-lb 


(17) 


TABLE VI 

OSCILLATOR INPUTS REQUIRED (VOLT) TO GENERATE PURE sin 44*- AND cos 41*- COMPONENTS 
OF PITCHING MOMENTS, ROLLING MOMENTS AND VERTICAL FORCES 


M 

^control * 

®os 

^oc 

^ss 

0$c 

^cs 

0 cc 

0.191 

^ control “ * 

+0.0143 


+0.0508 

+0.0290 


+0.0241 ; 


Me, control ~ * 

+0.0117 

lllfe 


+0.0283 


-0.0098 


^s, control “ 1 

-0.0177 

- 0.0236 

■ 

+0.0052 

-0.0169 

+0.0073 


L c, control ~ 1 

+0.0042 

-0.0071 

-0.0209 

-0.0200 

+0.0003 

- 0.0147 


T s , control ~ ^ 

+0.0922 

+0.1380 

-0.0490 

-0.0302 

+0.0252 

-0.0232 


T c, control ~ * 

-0.1044 

+0.1164 

+0.0123 

-0.0210 

+0.0235 

+0.0081 

0.239 

M s , control “ * 

+0.0028 

-0.0069 

+0.0299 

+0.0219 

-0.0111 

+0.0211 


M c, control “ * 

+0.0109 

+0.0028 

- 0.0096 

+0.0206 

■ 

-0.0070 


^s, control ” * 

-0.0023 

-0.0108 

-0.0056 

+0.0203 

I 

+0.0078 


L c, control = * 

+0.0128 

-0.0029 

-0.0245 

-0.0243 

-0.0008 

-0.0210 


\ control ~ * 

+0.1356 

+0.1337 

-0.0053 

-0.0155 

+0.0072 

-0.0128 


T c, control “ 1 

-0.1436 

+0.1085 

+0.0168 

+0.0070 

+0.0091 

+0.0100 

0.443 

M^ control ~ * 

- 0.0019 

rv rkrtc-o 
- U.UUJJ 

-i-n mcc 

+0.0! 16 

- 0.0069 

+0.0145 


Me, control = * 

+0.0053 

+0.0011 


+0.0331 

-0.0168 

-0.0004 


\ control “ * 

- 0.0057 

-0.0067 


+0.0253 

■ 

+0.0126 


L c. control “ * 

+0.0120 

- 0.0028 

■ 

-0.0084 

1 

-0.0112 


li 

i 

c 

8 

H* 

+0.1020 

+0.0732 

BeCo*’ iipl 

-0.0094 

-0.0018 

-0.0138 


T c, control * 

-0.0714 

+0.0941 

+0.0071 

-0.0108 

+0.0171 

- 0.0024 

0.849 

Ms, control = * 

+0.0049 

-0.0240 

+0.0338 

+0.0179 

- 0.0229 

+0.0182 


M c, control = * 

+0.0149 

-0.0149 

-0.0109 

+0.0487 

-0.0271 

- 0.0222 


^s, control = 1 

-0.0124 

-0.0137 

-0.0120 

+0.0074 

-0.0118 

+0.0024 


L c, control = * 

+0.0052 

-0.0056 

-0.0072 

-0.0121 

+0.0006 

-0.0123 


Ts, control — ^ 

+0.1050 

+0.0698 

-0.0211 

+0.0037 

+0.0017 

-0.0214 


Tc, control = 1 

- 0.0772 

+0.1079 

-0.0034 

-0.0305 

+0.0221 

+0.0031 

0.851 

n 

1 

c 

8 

Vi 

s 

+0.0001 

-0.0081 

+0.0191 

+0.0122 


+0.0109 


M c , control * 

40.0082 

- 0.0055 


+0.0290 


- 0.0055 


^s, control ~ * 

-0.0080 

-0.0107 


+0.0117 


+0.0098 


control “ 1 

+0.0113 

-0.0102 

- 0.0069 

+0.0028 


-0.0037 


Ts, control = * 

+0.1016 

+0.0599 

-0.0087 

+0.0109 


-0.0091 


Tc, control “ * 

-0.0682 

+0.0998 

+0.0034 

-0.0143 


- 0.0058 


* in.-lb 
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TABLE VII 

OSCILLATOR INPUTS REQUIRED (VOLT) 
TO COMPENSATE EXISTING 4P- VIBRATIONS 


TABLE VIII 

VIBRATION SUMMARY 


p 

0.191 

0.239 

0.443 

0.849 

0.851 

^os 

0.1683 

0.0394 

0.0146 

0.0457 

0.0300 

^oc 

0.3121 

0.0224 

-0.0490 

0.2354 

-0.2726 

^ss 

0.1746 

0.0090 

-0.1400 

-0.7980 

-0.3275 

^sc 

-0.0133 

-0.0293 

0.1273 

-0.5881 

0.3498 

^cs 

0.2052 

-0.0026 

-0.1 176 

0.4610 

-0.3549 

^cc 

-0.0651 

-0.0180 

0.0056 

-0.8308 

-0.0428 


Rel. Vibration 
Level 

M 

flnomi- 

nal 

Cj/o 

Ampl. of Pitch 
Variation 

1 

0.849 

10° 

-0.005 

^3.0° 

0.58 

0.851 

4 

-0.013 

2.0 

0.32 

0.191 

12 

0.102 

0.8 

0.21 

0.443 

4 

0.01 1 

0.5 

0.08 

0.239 

4 

0.028 

0.2 


Decreasing 

Vibration 

Level 


TABLE IX 

EFFECTS OF UNIT 4P OSCILLATOR INPUT ON BLADE BENDING 
AND TORSION MOMENTS (in-lb). p = 0.849 


M = 0.849 

Input 

sin 2^ 

cos 2^ 

sin 3^ 

cos 3^ 

sin 4^ 

cos 4^ 

sin 5^ 

cos 5^ 


^os 

0.3815 

- 2.6028 

- 1.1212 

+ 1 .9467 

+ 0.0022 

1.6252 

- 0.4640 

+ 0.2286 


$oc 

- 0.7265 

- 0.7428 

- 2.1170 

- 0.9082 

- 1.7646 

0.1744 

+ 0.4336 

- 0.2014 


^ss 

- 20.1796 

- 7.1252 

0.4843 

10.9746 

9.2290 

- 1.4705 

- 12.1221 

- 16.4408 

Flapbending 

^sc 

1.4455 

- 18.6069 

- 11.8793 

0.8771 

1.9670 

9.2946 

+ 18.4710 

- 13.1116 

3.3 in. 

^cs 

- 15.0717 

19.2091 

- 1.7568 

+ 13.3006 

4.4390 

13.0827 

24.2022 

- 18.4700 


^cc 

- 11.0041 

- 12.5052 

- 12.2451 

- 3.9250 

- 11.5818 

6.8481 

17.1269 

+ 18.2863 


0OS 

- 3.1446 

0.01156 

+ 0.0644 

- 6.4289 

0.5673 

- 5.5966 

- 2.6912 

- 5.2806 


0OC 

0.4488 

- 3.3139 

+ 5.7587 

- 0.6033 

7.2213 

1.7289 

4.4109 

- 1.9638 

Flapbending 

^ss 

- 13.1131 

- 1.6401 

- 9.4439 

11.4718 

2.7493 

1.6368 

20.3552 

30.4485 

13.15 in. 

0SC 

- 3.1093 

- 10.4663 

- 13.7168 

- 7.3647 

- 0.7250 

4.6008 

- 31.6355 

23.4534 


0CS 

- 15.3541 

3.9011 

- 20.8842 

- 14.1583 

- 4.0272 

- 4.6816 

- 53.1766 

36.9531 


fyx 

- 3.7738 

- 10.2279 

7.2742 

- 11.8491 

2.4534 

1.0036 

- 30.9619 

-33.3918 


&os 

- 5.2318 

5.1653 

18.4997 

- 66.4765 

8.5046 

- 2.0555 

6.0027 

8.6689 


0OC 

- 0.3311 

2.6008 

55.9170 

15.3823 

8.5503 

12.5308 

- 10.1401 

4.6381 


0ss 

- 23.2604 

3.6649 

-91.7693 

7.1537 

- 12.9172 

- 5.1116 

- 13.8450 

7.4174 

Chordbending 

^sc 

4.7043 

- 8.0015 

- 37.9514 

- 71.7419 

6.5301 

- 16.8130 

- 4.2184 

- 12.8505 

2.4 in. 

0cs 

- 25.0714 

15.3009 

- 59.7492 

- 177.5673 

41.5059 

-80.7110 

- 5.8153 

-27.4052 


^cc 

- 2.0059 

- 7.7253 

77.1483 

- 7.0902 

68.5358 

26.5134 

7.7451 

- 28.5566 


0OS 

0.1891 

0.0544 

- 0.2460 

0.5652 : 

- 1.0733 

0.2665 

0.1925 

0.0465 


doc 

0.0788 

- 0.1531 

- 0.1960 

- 0.2328 

- 0.6076 

- 1.0110 

0.0102 

0.01822 


^ss 

0.4975 

0.2685 

- 0.9271 

- 1.5838 

- 0.0498 

1 .4606 

15.6374 

13.1496 

Torsion 

^sc 

- 0.6976 

- 0.7498 

3.0700 

- 1.4345 

- 1.0039 

0.9952 

- 11.8807 

15.1709 

9.28 in. 

0CS 

0.8756 

- 0.0250 

- 1.5421 

- 0.9968 

- 1.9762 

1.0423 

- 14.6088 

21.3914 


dec 

- 0.8745 

- 0.9375 

2.1226 

- 2.5792 

- 1.3255 

- 1.2713 

- 17.9657 

- 13.8937 
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Comparison of Equations (15) and (17) shows that the vibratory 
pitching moment is reduced to approximately 15 percent and the 
rolling moment to approximately 18 percent of its original value. 
This indicates that the various blades behave differently and that 
the goal of zero 4P pitch-roll and vertical vibrations is achieved 
by cancellation of the effects of the four blades. 


Analytical Formulation and 
Calculated Results 

The aeromechanical characteristics of the High Advance 
Ratio Model (HARM) has been analytically described in 2 
degrees of freedom. These are based on the first and second flap- 
ping modes which have been approximated by polynomial fits of 


finite element determined mode shapes. The first and second 
mode shape approximations used are given by 

<t> x = 2.292x 2 - 1.292x 3 

and 

<t> 2 = -10.21x 2 + 20.78x 3 -9.57x 4 

where 

x = r^the non-dimensional radial station. 

R 


The aerodynamics are based on classical quasi-steady incom- 
pressible strip theory. The reverse flow region is fully accounted 
for, but stall effects have been neglected, as described m Refer- 
ence 5. 


TABLE X 

EFFECT OF VIBRATION COMPENSATION ON FLAPBENDING 
MOMENT (in-lb) AT 3.3 in. jj = 0.849 


n 


cos n 

sin n 

Amplitude 

2 

W/O Vibration Control 

-92.7652 

17.2338 

94.35 


Contribution of d Q 

- 0.0559 

- 0.1536 



0s 

16.6165 

15.2507 



0c 

19.2393 

2.2002 



TOTAL 

- 56.9653 

34.5311 

66.61 

3 

W/O Vibration Control 

- 1.1732 

- 14.7883 

14.83 


Contribution of 0 O 

- 0.1248 

- 0.5496 



fls 

- 9.2715 

6.5928 



6c 

9.3862 

9.3684 



TOTAL 

- 1.1833 

0.6233 

1.34 

4 

W/O Vibration Control 

- 0.1403 

- 3.5448 

3.55 


Contribution of 0 O 

0.1153 

- 0.4152 



0s 

- 4.2868 

- 8.5191 



0C 

0.3317 

11.6713 



TOTAL 

- 3.9801 

- 0.8078 

4.06 

5 

W/O Vibration Control 

3.2312 

2.2658 

3.95 


Contribution of 0 S 

- 0.0370 

0.0809 



0s 

20.8199 

- 1.1807 




-23.7042 

- 3.0926 



TOTAL 

0.3099 

- 1.9266 

1.95 
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TABLE XI 

SUMMARY OF OSCILLATORY BLADE LOADS (IN.-LB) 
WITHOUT AND WITH VIBRATION COMPENSATION 


Operating 

Condition 


M 

Flapbending at 3.3 

in. 

Flapbending at 13.15 in. 

Chordbending at 2.4 

in. 

Torsion at 9.28 in. 




n = 2 

n = 3 

u 

c 

n = 5 

n = 2 

n = 3 

n = 4 

n = 5 

n = 2 

n = 3 

n = 4 

n = 5 

n = 2 

n - 3 

n = 4 

n = 5 



0.191 

30.1 

4.4 

1.6 

3.5 

16.0 

1.9 

3.0 

4.3 

21.0 

2.2 

8.3 

19.4 

1.2 

0.7 

0.4 

0.6 

Without 


0.239 

10.5 

0.6 

0.2 

0.9 

5.3 

1.7 

0.9 

1.2 

4.6 

2.0 

11.0 

2.6 

0.5 

0.2 

0.3 

0.2 

Oscillatory ^ 


0.443 

16.4 

2.7 

0.1 

1.6 

9.2 

3.2 

0.4 

3.5 

9.4 

1.7 

10.5 

7.7 

0.9 

0.6 

0.3 

0.2 

Input 


0.849 

94.4 

14.8 

3.6 

4.0 

55.9 

3.6 

9.5 

5.9 

31.5 

31.4 

13.1 

14.6 

6.8 

4.1 

0.9 

0.3 


k 

0.851 

18.9 

8.6 

1.5 

3.1 

17.7 

4.6 

3.4 

5.8 

17.4 

10.9 

18.9 

10.7 

3.3 

2.4 

0.7 

0.4 


r 

0.191 

29.6 

1.1 

2.9 

0.4 

16.1 

4.4 

5.0 

3.0 

19.2 

22.7 

10.9 

3.9 

1.0 

1.2 

0.4 

4.5 

With 


0.239 

10.3 

0.4 

0.3 

0.7 

5.3 

1.9 

0.8 

1.5 

4.7 

3.0 

11.5 

2.1 

0.5 

0.3 

0.3 

1.2 

Oscillatory 

Control 


0.443 

12.3 

1.3 

1.3 

1.1 

7.5 

2.7 

0.5 

1.3 

7.7 

3.5 

13.6 

8.7 

0.8 

1.4 

1.4 

1.7 

Input 


0.849 

66.6 

1.3 

4.1 

2.0 

41.7 

1.3 

2.4 

2.1 

15.7 

68.8 

38.9 

22.3 

6.5 

4.4 

0.8 

3.1 


k 

0.851 

20.2 

2.4 

6.5 

4.1 

16.5 

3.8 

7.0 

2.5 

17.5 

13.6 

75.6 

7.0 

2.3 

4.0 

0.7 

3.5 


The method of solution provides for intermode harmonic 
coupling through the 17th harmonic. This is accomplished by 
obtaining transient solutions of the 2-degree-of-freedom descrip- 
tion of the rotor system described as constant coefficient linear 
differential equations over 10-degree sectors of the rotor 
azimuth. 

The values of the coefficients for the system of differential 
equations evaluated in this work have been determined at the 
center of the sectors i. e., at 5°, 15°, 25°, etc. 

The basis for the analytical formulation is founded on 
Shannon’s sampling theorem which says that the discrete signal is 
equivalent to the continuous signal, provided that all frequency 
components of the latter are less than 1/2T cycles per second, T 
being the time between instants at which the signal is defined, 
(References 6 and 7 ). Since the solution also provides for a com- 
pletely general transient solution, it can be used to calculate a 
Floquet solution by specializing the initial conditions. This has 1 
been done for the square spring oscillator case studied by M. A. 
Gockel and reported in the AHS Journal in January 1972. The 
problem statement which is exactly describable by this theo- 
retical method was shown to yield the identical Floquet solutions 
as those reported. It is important to note that should the system 
be unstable, the harmonic balance method of solution would not 
directly reveal this instability. 

Briefly, the initial conditions at the beginning of a sector are 
determined by calculating the terminal conditions for the pre- 
vious sector which are then used to initialize the new sector. It 
has been found that essentially arbitrary conditions can be used 
to start the solution and that excellent steady-state conditions 
have been obtained for the conditions examined in six rotor revo- 
lutions. For each solution case presented, the rotor has been 
solved for eight revolutions to ensure that the second flapping 
mode contribution to the response has converged to a steady- 


state value accurate to at least four significant figures. The pro- 
gram is used to calculate closed-form analytic solutions over each 
10-degree sector and therefore is not dependent on a particular 
method of numerical integration. (See Appendix A.) The 
method, however, when applied to the analysis of steady-state 
conditions, does require that sufficient solution time be calcu- 
lated so that initial transients are dissippated to ensure that 
steady-state equilibrium is achieved (Reference 8). 

The test configuration experimentally examined with re- 
spect to IP flapbending distributions at p = 0, including center- 
line measurements, has been compared with this analysis 
procedure on Figure 2, utilizing the two-mode description. This 
is a limited use of the analysis technique to establish test/analysis 
correlation. It is believed that the absence of time-dependent 
aerodynamics quasi-steady, largely accounts for the phase error 
in response. The centerline shaft moment measured was 0.75 of 
the calculated (a = 5.73). This may be due to the relatively low 
inflow of the test condition. 

In general this correlation, including the spanwise distribu- 
tion, appears reasonable. 

The eigenvalues of each 10-degree sector are evaluated as 
part of the method. These are summarized in Tables XII, XIII, and 
XIV versus azimuth the /i = 0. 191 , 0.45, and 0.85 where the real 
and imaginary parts of the eigenvalues have been normalized by 
the noted natural-mode frequencies. The negative aerodynamic 
spring effects over azimuth 90 < V F< 270 as well as the positive 
stiffening from 270 < <90 are as expected more pronounced 

on the first mode frequency. The effects of reduced aerodynamic 
spring and damping are also seen on the retreating side. These 
results show that both damping as well as frequency variations 
occur around the azimuth which influence the rotor response 
with harmonic excitations. 
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3.0 


TABLE XII 

NORMALIZED EIGENVALUES* AT EACH 10-DEGREE 
AZIMUTHAL SECTOR FOR p = 0. 1 9 1 




Figure 2. One-Per-Rev Blade Radial Flap-Bending Moment 
Distribution at p- 0. 


The rotating frequencies and properties of the flapping 
modes noted in Tables XII, XIII, and XIV analytically describe 
the 7.5-ft-diameter rotor, configuration (5), 500-rotor-rpm con- 
dition for which all harmonic feathering tests were conducted. 
In an effort to further improve analytic correspondence with 
test data the slight change of the second flapping mode fre- 
quency resulted from matching collective blade angle selection 
at the test conditions. Details of the test model are given in 
References 9, 10 and 1 1 . 



♦SECTOR EIGENVALUES ARE GIVEN BY: 


The harmonic components of the blade root flap-bending 
moment (OP through 5P) were calculated for these advance ratios 
for unit perturbation of blade feathering angle at 0 lc , 0 Js , 0 2c , 

#2s' ^3c’ e 3s’ ^4c’ ^4s’ ^5c’ ^Ss 1 as wel1 as for un * t chan 8 e * n 
0 O and a 


The single non-dimensional blade root, centerline flap-bending 
moment harmonic influence coefficients resulting from harmonic 
feathering are summarized in matrix form in Tables XV, XVI, and 
XVII for n = 0. 1 91 , 0.45, and 0.85. These are based on har- 
monic analysis of the moment at each condition for 36 equally 
spaced ( 1 0-degrees apart) azimuth intervals. Single-blade 


(Rj +l x i) ( 1.3412 ) 

AND (R 2 + I 2 i) (6.3812 ) 

computed root flap-bending moment influence coefficients 
at H = 0.45 are compared with experimental 0.073R 
single-blade data, in parentheses, from Reference 1 and 1 3 
in Table XVIII. 

These appear reasonable when shear effects are considered. 

It is important that the general character of these influence 
coefficients be established in future tests. These tests should be 
structured to permit measurement to confirm these distributions. 
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♦SECTOR EIGENVALUES ARE GIVEN BY: 
(Rj+1,0 (1.3412 ) 
AND (R 2 + I 2 0 (6.20 0 ) 


♦SECTOR EIGENVALUES ARE GIVEN BY 
(Rj ± i! i) (1.34 fi ) 
AND (R 2 ± 1 2 J2) (6.20 n) 



c TABLE XV 

-ML- BLADE ROOT (STA 0) BENDING MOMENT INFLUENCE COEFFICIENT MATRIX FOR v = 0 191 

(Pj = 1.34, P 2 = 638) 




CPo 

c Pic 

CPis 

Act 

.0034 

.0009 

-.0016 

A0 

.0132 

.0049 

-.0111 

A0JS 

.0036 

.0057 

-.0225 

A01C 

-.0005 

-.0213 

-.0045 

A02S 

0 

-.0056 

-.0004 

A02C 

-.0003 

-.0005 

.0057 

A03S 

0 

-.0001 

.0004 

A03C 

0 

.0004 

0 

A94S 

0 

0 

0 

A64C 

0 

0 

0 

A05S 

.0119 

0 

0 

A05C 

0 

0 

0 


CPjc 

C P2S 

C P3C 

C P3S 

-.0001 

-.0001 

0 

0 

-.0007 

-.0005 

-.0001 

0 

-.0018 

-.0003 

0 

.0002 

-.0001 

.0018 

.0002 

0 

.0018 

.0031 

.0014 

-.0009 

.0031 

-.0018 

-.0009 

-.0014 

.0001 

.0002 

-.0046 

-.0012 

.0002 

-.0001 

-.0012 

.0046 

0 

.0001 

-.0010 

.0017 

.0001 

-.0001 

.0017 

.0009 

0 

0 

.0002 

.0002 

0 

0 

.0002 

-.0002 


C P4C 

0 

0 

0 

0 

.0002 

.0003 

.0015 

-.0014 

-.0076 

-.0026 

-.0015 

.0024 


C P4S 

0 

0 

0 

0 

.0003 

-.0002 

-.0014 

-.0015 

-.0026 

.0076 

.0024 

.0015 


C P5C 

0 

0 

0 

0 

0 

0 

.0001 

.0003 

.0017 

-.0020 

-.0101 

-.0033 


C ft5S 

0 

0 

0 

0 

0 

0 

.0003 

-.0001 

-.0020 

-.0017 

-.0033 

.0101 


r TABLE XVI 

-Mi- BLADE ROOT (STA 0) BENDING MOMENT INFLUENCE COEFFICIENT MATRIX FOR M =.45 

(Pj = 1.34, P 2 = 6.20) 



cp 0 

epic 

CPlS 

C P2C 

c p 2S 

C P3C 

C P3S 

C P*1C 

C P4S 

c Psc 

c Pss 

A a 

.0085 

.0053 

-.0089 

-.0010 

-.0012 


-.0004 

-.0001 

-.0001 


0 

A0 

.0160 

.0135 


-.0038 

-.0029 


-.0004 

-.0001 

-.0002 


0 

A01S 

.0087 

.0102 


-.0048 

-.0016 


.0007 

-.0006 

0 


0 

A01C 

-.0011 

-.0226 

-.0034 

-.0004 

.0046 

.0011 

.0002 

.0001 

.0003 


0 

^02S 

0 

-.0130 

-.0006 

.0006 

.0046 

.0036 

-.0020 

.0009 

.0015 


.0002 

A02C 

-.0015 

-.0024 

.0142 

.0043 


-.0020 

-.0035 


-.0009 


.0002 

A03S 

0 

0 

.0023 

.0004 

.0010 

-.0057 

-.0032 

.0038 

-.0036 

.0008 

.0016 

A03C 

-.0002 

.0022 

-.0002 

.0007 

-.0001 

-.0033 

.0057 

-.0036 

-.0038 

.0016 

-.0008 

A04S 

0 

-.0001 

0 

.0003 

.0007 

-.0026 

.0042 


-.0046 

.0043 

-.0048 

A64C 

0 

0 

0 

.0004 

0 

.0042 

.0026 

-.0046 

.0090 

-.0048 

-.0043 

A05S 

0 

0 

0 

0 

.0003 

.0008 

.0008 

-.0038 

.0058 

-.0118 

-.0055 

A65C 

0 

0 

0 

0 

.0003 

.0008 



.0038 

-.0054 

.0118 


C RM TABLE XVII 

- BLADE ROOT (STA 0) BENDING MOMENT INFLUENCE COEFFICIENT MATRIX FOR fi = 85 

(Pj = 1.34, P 2 = 6.20) 


A a 

A6 

A61S 

A01C 

A02S 

A02C 

A03S 

A03C 

A04S 

A04C 

A05S 

A05C 




'-Pic 


KjMjf 

.0227 


.0378 


.0278 


-.0258 

i 

-.0229 


-.0110 


-.0014 

| 

.0060 

i 

-.0004 

-.0004 

.0002 

0 

.0001 

0 

.0006 


'-Pis 


^P2C 


W*2S 


^P3C 


L P3S 


C P 4c 


C P4S 


C P5C 


c h$ 


-.0296 


-.0102 

-.0598 

-.0141 

-.0155 


-.0117 

-.0114 


-.0012 

.0085 


-.0026 

.0079 


.0084 

.0067 



.0035 

0 

.0029 

-.0009 


.0013 

.0008 


.0007 

-.0014 



-.0003 


-.0003 

-.0002 


-.0037 

-.0061 

-.0036 

.0035 

.0081 

-.0026 

-.0088 

-.0084 

-.0067 

.0087 

.0029 

.0023 


-.0039 

-.0039 

-.0004 

.0008 

-.0034 

-.0064 

-.0082 

.0089 

.0087 

.0067 

.0023 

-.0029 


0 

0 

-.0019 

.0003 

.0024 

.0052 

.0100 

-.0084 

-.0135 

-.0100 

-.0088 

.0124 


-.0021 

-.0032 

-.0014 

.0022 

.0054 

-.0019 

-.0084 

-.0100 

-.0100 

.0134 

.0124 


0 

-.0003 

-.0005 

0 

-.0019 

.0013 

.0033 

.0048 

.0112 

-.0101 

-.0167 


-.0015 

-.0018 

-.0001 

.0009 

.0016 

.0023 

.0049 

-.0033 

-.0100 

-.0112 

-.0115 
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TABLE XVIII 

BLADE ROOT (STA 0) BENDING MOMENT (IN-LB)/DEG INFLUENCE MATRIX FOR ji = 45 

( 52.36, Pj - 1.34, P 2 = 6.20) 



Po 

Pic 

Pis 

P2C 

P2S 

P3C 

P3S 

P4C 

P4S 

P5C 

Pss 

LIFT 

Aa 

19 

12 

-20 

-2 

-3 

-1(1) 

-1 (-2) 

0 

0 

0(1) 

0 (-1) 

6 

A0 

36 

31 

-62 

-9 

-7 

-2(1) 

-1(1) 

0 

0 

0(1) 

0(0) 

10 

A01S 

20 

23 

-66 

-11 

-4 

-1 (1) 

2 (-2) 

-1 

0 

0(1) 

0(0) 

6 

A81C 

-2 

-51 

-8 

-1 

10 

3(0) 

0(1) 

0 

1 

0(0) 

0(0) 

0 

A02S 

0 

-29 

-1 

1 

10 

8 

-5 

2 

3 

-1 

0 

0 

A02C 

-3 

-5 

32 

10 

0 

-5 

-8 

3 

-2 

0 

0 

-1 

A03S 

0 

0 

5 

1 

2 

-13 

-7 

9 

-8 

2 

4 

0 

A03C 

0 

5 

-1 

2 

0 

-7 

13 

-8 

-9 

4 

-2 

0 

A04S 

0 

0 

0 

1 

2 

-6 (0) 

9(6) 

-20 (-8) 

-10 (-5) 

10 (-5) 

-11 (-1) 

0 

A04C 

0 

0 

0 

1 

0 

9(6) 

6 (-2) 

10 (-4) 

20 (7) 

-1 1 (-6) 

-10(3) 

0 

A05S 

q 

0 

0 

0 

1 


O 

A 

-9 

13 

~is / 

-12 

A 

KJ 

A05C 

0 

0 

0 

0 

1 

2 

-2 

13 

9 

-12 
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Full-Scale Control Loads 

The feasibility of active vibration attenuation depends on 
the capability of the rotor to generate cancelling shaft moments 
and shears while control forces and displacements remain within 
acceptable limits. 

Since full-scale data are the most relevant from the stand- 
point of hardware test background, the CL 840/AMCS 
(Advanced Mechanical Control System) Cheyenne rotor 
configuration, at a gross weight of 20,000 and with a rotor shaft 


moment of 100,000 in.-lb, was analyzed for hovering flight to 
gain a numerical measure of how loads compare with limits. In 
this analysis three higher harmonic blade-feathering excitations, 
3P, 4P and 5P, were examined to determine the relationships 
among control loads, shaft moments and shear forces. The 
Lockheed Rotor Blade Loads Prediction Model was used for 
this analysis; 68 finite elements were used to describe the system. 
The calculated results, based on 1 -degree excitation levels, are 
summarized in Table XIX. 


TABLE XIX 
CL 840 ANALYSIS - 

SHAFT AND BLADE LOADS DUE TO ONE-DEGREE 
OF HIGHER HARMONIC BLADE-FEATHERING MOTIONS 





FEATHERING FREQUENCY 




3^ 

4 <P 

5* 


Endurance 


Amplitude 

Phase 

Amplitude 

Phase 

Amplitude 

Phase 

Limit, in.-lb 

Shaft Forces 








4P H-force 

380 lb 

61° 

40 1b 

59° 

310 1b 

34° 


4P Y-force 

380 lb 

84° 

40 lb 

83° 

3101b 

12° 


4P Pitching Moment 

22,000 in.-lb 

83° 

0 


108,000 in.-lb 

8° 

| 325,000 

4P Rolling Moment 

22,000 in.-lb 

16° 

0 


108,000 in.-lb 

76° 

4P Thrust 

0 


3000 lb 

o 

o 

0 



Blade Root Torsion * 








Harmonic 








Steady 

-3800 in.-lb 


-4000 in.-lb 


-3900 in.-lb 



IP 

210 in.-lb 

11° 

210 in.-lb 

11° 

220 in.-lb 

11° 

) 

2P 

80 in.-lb 

49° 

50 in.-lb 

42° 

50 in.-lb 

39° 

l 15,500 

3P 

1500 in.-lb 

15° 

70 in.-lb 

82° 

40 in.-lb 

OO 

A 

O 

4P 

130 in.-lb 

47° 

13300 in.-lb 

O 

00 

00 

400 in.-lb 

57° 


5P 

20 in.-lb 

27° 

80 in.-lb 

35° 

7800 in.-lb 

10° 



* Pitch link forces are internal loads between the blade and 
swashplate and therefore sclf-cancelling. 


274 



The calculated root torsion moments shown in the table 
reflect both the feathering moments at the primary exciting 
frequency and the interharmonic coupling terms; as expected, 
the latter are considerably less. Pitch link loads can be 
determined by multiplying the root torsion moment by 0. 1 (to 
account for all applicable geometry): endurance limit of the 
pitch link load is 1550 pounds. 


The 7.5-foot hingeless rotor model data showed that 0.2 to 

0. 6-degree cyclic angle excitation levels were required. Study of 
CL 840 test data indicates that similar blade excitation would 
be expected with a full-scale, four-bladed rotor. The CL 840 
data are not yet published in documents that can be referenced, 
however, this material is expected to be published during 1974. 


In summary, full-scale data founded on endurance limit 
considerations indicate that internal blade loads and control 
loads will not limit the trim flight use of periodic variation of 
conventional controls for vibration attenuation. 


Conclusions 


The present report is a preliminary evaluation of the con- 
cept of vibration reduction by properly selected oscillatory col- 
lective and cyclic control applications. The investigations are 
based on experimental frequency response data covering advance 
ratios from approximately 0.2 to 0.85. 

Because there was no instrumentation for the measurement 
of the pitch and roll vibrations, these values were obtained by 
properly adding up the flap-bending moments at 3.3 inches. Any 
other quantity representing pitch/roll vibrations can be 
compensated for in the same fashion. 

The calculated control inputs required for vibration reduc- 
tion stay within acceptable limits. For four of the five conditions 
tested they are smaller than the values used for the frequency 
response tests. The blade pitch variations required for vibration 
alleviation vary, depending on the advance ratio, less than 1° for 
.2 p < .45 and - 3° for p = .85. 

As to be expected, the compensating controls greatly affect 
the blade loads, Le., torsion, flap- and chordwise bending. With 
regard to flap-bending at 3.3 inches (root flexure), the following 
statements can be made: 


• 3 and 5P flap moments were, by command, 
drastically reduced 

• 2P flap moments were least affected. These 
were the largest oscillatory loads. 

• 4P flap moment increments generally increased 
with increasing advance ratio, but were small 
relative to the 2P flap moments. 

As a general rule, chordwise bending and blade torsion 
increments also increase with the advance ratio. At lower (i 
values the loads are not critical. It is concluded that the 
concept investigated is primarily suited for low and medium 
advance ratios, i.e., for the speed-range of present day 
rotary wing aircraft. The latter application appears promising 
and further studies and tests are suggested. Instrumentation 


to determine rotor vertical and inplane shear forces should 
be incorporated in such future tests. Also a system with a 
first inplane frequency in the vicinity of 1 .5P in combination 
with a flapping frequency of 1.1P should be tested at con- 
ventional advance ratios to provide experimental data 
representative of current designs. 
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where o^, o 2 , • ■ ■ • <* n are all distinct, this yields 


(3(t) = 


n 

2 

k=l 


r(V 

Q>k> 


e 


“kt 


In the case cited 


12. Watts, G. A. and London, R. J., VIBRATION AND LOADS Q( s ) = A(s) (s - a) (s - y) 

IN HINGELESS ROTORS, Vol. I and II, NASA 

CR-1 14562, September 1972. where = 0 

a \ u 


( 4 ) 


Appendix A 

The transient response solution of a system described by 
constant coefficient linear differential equations is developed 
in this appendix. The single-degree-of-freedom case with arbi- 
trary initial conditions and solution of the general case for an 
nth order system with both zero and nonzero initial conditions 
is reported. 


a 2 ~ ot 
o 3 = y 

P(s) = L + [ P(0)A] s 2 + [ p(0)B + P(0)A] s 
Therefore 


Given the single degree of freedom: 

A^4+ B^f + CP = F(t) 
dt 2 dt 

where A, B, and C are constants, then 

Az(^f) +BjC^)+CX<P) = X(F(t)) 
where Ji is the Laplace transform operator. This yields 


p(t) “ k ?, W e 


pa) = 


A(- »)(- y) 


L + [p(0)A] o 2 + [p(0)B + p(0)A] o 


A(+ o') (+ a - y) 


c 0 ‘ (5) 


L + [(3(0) A] Y 2 + fp(0)B + p(0)Al y 


(A)(+y)(y- a) 


(As 2 + Bs + C) p(s) = F(s) + p(0)(As + B) + P(0)A (2) 

or 

s) = F(s) + p(0) (As + B) + P(0)A 
As~ + Bs + C 

If a positive constant step load of magnitude + L is the form of 
F(t), then 


Extension to the general case is accomplished as follows. 
Given the general determinantal equation: 

{s 2 ( a] + s [ B ] + [c]} jp(s)j = |f(s)| (6) 

Where the elements of matrix A, B, and C are constants, 
using Cramer’s Rule: 


jC(F(t)) = F(s) = ~ 


and 


- L p(0)As 

A(s)(s - o-Hs - y) A(s-a)(s-y) ^ Expanding 

+ P(0)B + p(0)A 

Afs-oOCs-y) yields 

Where p(0) and j3(0) are the values of the variable pat time 
t = 0 and a, y are the roots of s^ + Bs/A + C/A, p(s) trans- where 

formed back into the time plane is accomplished through use 
of the inverse Laplace transform of the form P(s) 

Q(s) 


where 

P(s) = polynomial of degree less than n 


Denominator with 

Column i replaced by F(s) (7) 



A 0 (s - o 1 )(s - a 2 ) .. • (s-cr n ) (8) 

A q = Coefficient of highest power term 

a- (i = 1 . . . n) are the eigen values (roots) of 
the determinantal equation 


and 


Q(s) = (s - a j ) (s - a 2 ) . . . . (s - a n ) 
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Case 1 - Zero Initial Conditions 

Assume p-(0) and p-(0) for ail i are both zero and that a 
positive unit load acts on p e and that the response of pf is to 
be determined. Then 


Case 2 - Nonzero Initial Conditions 

The general form of F(s) now becomes: 


F(s) } = { + 1 /s in row e with all other 

rows equal _0 


F(s) = 


+ s A + 


h |[A]jp i( 0)jJ 


[B]Jj Pi (0) 


(ID 


Defining 


A + s B + C 


(e,f) 


(9) 


as the original determinantal equation with Row e and Column 
f removed and all the remaining rows and columns moved up 
and to the left, respectively, this forms a determinantal equation 
of one less order. 


where Lj are the forces applied at each coordinate p- and p-(0) 
and j3j(0) are the positions and rates of the coordinates at time 
zero (initiations of the solution). In this case place the column 
s M into the column location of the coordinate for which the 
response is desired without reduction of the order. Then 


Based on the earlier development in the s-plane 


Column i 



X. M) =^ + JL + ^L + ... + ^i_ 

^ V 1 f S S-ttj s-o^ s -<* n 

P(s) = 

s|F(s)J 


(12) 

and in the time plane 

where all other terms are 



„ , x a lt oh t at 

3f(t) - a Q + a j e 1 + a 2 e 1 + . . . + a n e n 

< 10 > | s 2 [ A J 

+ S ( B I + ( C I 1 


where 

(_l)(e + D D(o) e f 
3 ° n 

and 




Q(S) = 

Aq(S-Oq) (s-g 


(13) 


and 


A 0 n^j 
i=l 

_ (.pCe+QD^ f 
n 

»j a q n -<>j) 

i=l 


i t j 


wjiere the a’s are the eigenvalues of the determinantal 
equation 


s|s 2 [a] + s[ B ] + [c] |J - 0 


A Q is determined by the relationship 


D (o) = a 0 n «i 

i=l 


Then 


v P(o^) t> k t 


(14) 


D(o) e f an< * D(a:) e f are f Qrine d from the original determinantal 
equation with Row e and Column f removed and all the remaining 
rows and columns moved up and to the left, respectively, 
evaluated at o andcxj. The £*j are the roots of the original deter- 
minantal equation before Row e and Column f were removed. 
These roots are assumed distinct, an unimportant limitation for 
most physical systems. Note that this solution does not preclude 
instability either aperiodic or oscillatory. 

In practice the eigenvalues are obtained prior to the forma- 
tion of the coefficients and are examined to verify the distinct 
character of the eigenvalues. 

Scalar multiplication of this solution provides the result 
for the nonunit loading case. Summation of solutions obtained 
for loadings at each coordinate can be used to provide the 
general solution for this case where p^O) and ^(0) for all i 
are both zero, i.e., that the initial conditions at time zero 
are all zero. 


where s “ 0 and the remaining eigenvalues of the general deter- 
minantal equation form the set of a k ’s, and Aq is determined 
by the relationship 

IHO) = A 0 naj (15) 

i— 1 

as given in Case 1 . 


In most applications the restriction that the initial condi- 
tions are zero is an unacceptable constraint and this condition 
has been relaxed; the solution follows. 
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